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ABSTRACT 



We consider for the first time the implications on the modified gravity MOND model of galaxies, of the presence of 
dark baryons, under the form of cold molecular gas in galaxy discs. We show that MOND models of rotation curves 
are still valid and universal, but the critical acceleration ao separating the Newtonian and MONDian regimes has a 
lower value. We quantify this modification, as a function of the scale factor c between the total gas of the galaxy and 
the measured atomic gas. The main analysis concerns 43 resolved rotation curves and allows us to find the best pair 
(ao = 0.96 X 10~^° m.s~'^, c = 3), which is also compatible to the one obtained from a second method by minimizing 
the scatter in the baryonic TuUy-Fisher relation. 
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1. Introduction 

The missing mass problem at galactic scales is revealed by 
rotation curves of spiral galaxies, where stars and gas ro- 
tate with speeds higher than expected (e.g. Sofue & Rubin, 
2001), in the frame of Newtonian gravity. This missing 
mass has been called dark matter and has now been de- 
veloped into a standard cosmological model ACDM (e.g. 
Blumenthal et al. 1984). Dark matter particles collapsed 
to form the first structures of the Universe, and it is in 
those potential wells that the baryonic gas infalls and cools 
to form galaxies. As a result, galaxies are embedded in a 
spheroidal dark matter halo with a mass profile obeying 
M(r) cx r in the outer parts, so that the Keplerian circu- 
lar velocity, Vc = \/GM[r)lr^ tends to be constant in the 
majority of observed rotation curves (e.g. Bosma 1981). In 
the central parts however, the CDM numerical simulations 
predict a cuspy dark matter profile, with a relative CDM 
amount much larger than what is observed (Navarro et al 
1997, Navarro & Steinmetz 2000). 

An alternative, without invoking any dark matter, is 
to consider a modification of the Newtonian gravity law 
at low accelerations, so that the luminous mass is suffi- 
cient to describe the dynamics of the gravitational system. 
This was done empirically by Moti Milgrom who proposed 
the Modified Newtonian Dynamics (MOND) paradigm 
(Milgrom, 1983). The MOND regime is distinguished from 
the Newtonian regime when the characteristic acceleration 
of the gravitational system falls below a critical accelera- 
tion ao 1-2 X 10"^" m.s~^. In the deep MOND regime, the 
modification au of the Newtonian acceleration, oat, can be 
written as om = y/cioo-N- The asymptotic circular velocity 
is thus a constant : v"^ — GMqq. This is the TuUy-Fisher 
law, which finds there a justification. 

In this paper, we consider for the first time how the 
problem of missing baryons could be made compatible with 
the MOND phenomenology. In the cosmic baryon budget. 



it is now well known that only about 6% of them shine as 
stars and gas in galaxies (Fukugita et al. 1998). Locally, 
UV absorption in front of background sources have shown 
that as much as 30% of baryons could be associated with 
the Lymana forest and an uncertain fraction (5-10%) with 
the warm-hot medium (WHIM at lO^-lO'^K, Nicastro et al 
2005, Danforth et al 2006). So about half of the baryons are 
not yet accounted for. Most of them should be in the cos- 
mic filaments, in the intergalactic medium, but it appears 
unavoidable that a significant fraction of them exist as cold 
gas in galaxies themselves, as suggested by Pfenniger et al 
(1994), and Pfenniger & Combes (1994). From the rotation 
curves, it is possible to put an upper limit on the frac- 
tion of dark baryons in galaxies: the dark baryons cannot 
overpredict the rotation speeds. HI interferometric stud- 
ies of nearby galaxies have shown that the dark matter 
distribution necessary to fit the rotation curves of galax- 
ies follows the surface density of the atomic gas (Bosma 
1981). Hoekstra et al. (2001) estimate that a factor of 7-10 
between the dark baryons and atomic gas reproduces the 
rotation curve of any galaxy, giant or dwarf, high or low sur- 
face brightness, without adding any dark matter halo. This 
scale factor means than there cannot be more than 10% of 
the total baryons in the form of cold gas in galaxies. The 
observed HI extended rotation curves roughly require only 
to double the fraction of baryons in galaxies at maximum. 
Out of the 50% of the dark baryons still not accounted for, 
no more than 10% can be in galaxies, but they can be less 
abundant. 

Another way to trace the dark baryons in galaxy discs 
is the TuUy-Fisher relation. This observational law involves 
the luminosity versus rotational velocity, which is equiv- 
alent to the stellar mass versus the velocity (with a con- 
stant mass-to- light ratio). A divergence from this law was 
found for late-type galaxies. On the TuUy-Fisher diagram, 
the stellar mass of these galaxies is too small for the ve- 
locity observed. The problem was solved by taking into ac- 
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count the mass of the atomic gas instead of only the steUar 
mass (McGaugh et al., 2000). The next step was to sup- 
pose that all the gas is not seen (dark baryons). By this 
approach, Pfenniger & Revaz (2005) find that a factor of 
3 between the total and atomic gas reduces the scatter of 
the baryonic TuUy-Fisher relation. More recently, Begum 
et al. (2008) carried out the same analysis but focusing on 
a sample with low-mass galaxies, which are more constrain- 
ing because they generally contain more gas. They find a 
scale- factor of around 9. 

From a dynamical point a view, Revaz et al. (2008) per- 
formed some numerical simulations of galaxies in Newton 
gravity with dark matter. They add a cold dark baryon 
component and show that a factor of up to 5 between the 
total and atomic gas is realistic to reproduce the global 
behavior of the galaxy: the stability against local axisym- 
metric collapse is ensured through a low dissipation, and 
patterns still spiral can develop due to the self-gravity of 
the disc. 

It is thus interesting to take into account the possible 
existence of dark baryons in the framework of MOND at 
galactic scales. If some fraction of the assumed dark mat- 
ter is real, the actual critical acceleration «□ is overesti- 
mated. We want to derive the new possible limiting accel- 
eration, and interpolating function fi. The contribution of 
dark baryons has been discussed in modeling galaxy cluster 
cores by Milgrom (2007). At galactic scales, Gentile et al. 
(2007) have studied the contribution of neutrinos to rota- 
tion curves. 

The paper is organized as follows: in section 2, we in- 
troduce the contribution of dark baryons at galactic scale, 
in section 3 we describe the galaxy sample. The modeling 
of galaxies is discussed in section 4. In section 5, the impli- 
cations of the existence of dark baryons with MOND are 
explored through rotation curve fits and the TuUy-Fisher 
law, then discussed in section 6. 



2. Overview 

We consider that part of the missing mass problem in 
galaxies could be alleviated by the presence of some dark 
baryons. They could be in the form of cold molecular gas 
II2 in the outer flaring disc. This scenario has been stud- 
ied in the framework of Newtonian gravity (Pfenniger & 
Combes, 1994; Combes & Pfenniger, 1997). In this model 
the amount of non-baryonic DM condensed in galaxies is 
reduced, but still necessary. 

At present, in the framework of modified gravity, since 
the dark molecular gas increases the baryonic mass, the 
critical acceleration of MOND will be reached at a lower 
level, and this is not necessarily in agreement with a sin- 
gle acceleration qq for all galaxies. Is the presence of some 
dark baryons in galaxies still compatible with MOND? As 
the required "phantom" dark matter is proportional to the 
HI gas (Hoekstra et al 2001), it appears quite possible to 
combine the presence of dark baryons with MOND, by only 
reducing the value of the critical acceleration. For example, 
we examine the dwarf LSB galaxy NGC 1560, dominated 
by dark matter in the Newtonian frame. Its rotation curve 
reaches a plateau at 7 kpc. In accordance with MOND, 
this maximum velocity can be written: 

v'^ = GMao. 




Radius (kpc) 

Fig. 1. Rotation curve fit of NGC 1560, in the MOND 
framework with dark baryons in the form of cold molecular 
gas. If this molecular gas is such that mass Mgas = cMat, 
the value of the actual estimation of the critical accelera- 
tion ao = 1.2 X 10^^" m.s~^ must be reduced (divided by a 
factor c, to a first approximation) , to be in agreement with 
the observations. The factor c is indicated in the bottom 
right corner, for each curve fitted. 

The mass M, to a first approximation, corresponds to the 
mass of the atomic gas (HI-|-He) Mat, since at this outer 
radius, the stellar content is negligible. If we consider that 
this galaxy may contain as much molecular gas as atomic 
gas, which means that the total mass is Mgas = '^Mat, the 
actual value oq must be divided by 2 to conserve the same 
asymptotic velocity. With this model, if the total mass of 
gas is equal to the atomic gas mass times a scale factor, 

Mgas = cMat, 

the critical acceleration uq must be divided by c to keep the 
same asymptotic velocity. It works correctely for the galaxy 
NGC1560, as shown in the Fig. [T] The most remarkable 
is the shape of the rotation curve which is conserved by 
modifying the values (c, ao), not just the maximum velocity. 
The wiggle at 5 kpc is well reproduced until c = 5. Beyond 
this value, the oscillation is more and more pronounced and 
no longer matches the data. 

It might be surprising that two very different models 
succeed in representing equally well the rotation curves of 
galaxies: the first is to assume the presence of 7-10 times 
more gas than observed in HI, the second is to modify the 
gravity. For dwarf irregular galaxies, where the stellar com- 
ponent is minor, the visible matter is dominated by the HI 
gas and the total mass is dominated by the dark matter. 
With these approximations, we can write two equations: 
in the Newton model, the acceleration at the outskirt of 
the disk is V^^^/r ^ 7GM{HI)/r'^, since the total mass is 
very close to 7 times the gas mass for these dwarfs, and 
in the MOND model, V^^Jr - {aoGM{HI)/r^y/^ , since 
they are in the deep MOND regime in the outer parts. 
Combining these two relations implies that the accelera- 
tion in these dwarfs is always about 7 times less than ao, 
which is remarkable (cf also Fig [3]). As for the radial distri- 
bution, even for the dwarf galaxies, the regime is still par- 
tially Newtonian in the center, and the two models work 
equally well, because of a combination of factors including 
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morphological type ^, (kms^^) 

Fig. 2. Left: Morphological type in the sample made of 43 
nearby galaxies. Right: Distribution of the maximum veloc- 
ity in the galaxy rotation curves. 

the stellar component and the interpolation MOND func- 
tion (see Figure 1). 

In this paper, we study whether MOND is still consis- 
tent with the presence of dark baryons such as the molecu- 
lar gas H2. Is there still a critical acceleration universal for 
all the galaxies? What is the proportion of dark baryons 
that is consistent with the observations, in the framework 
of MOND? We try to answer these questions in two ways. 
First we analyze a first sample of galaxies by fitting their 
rotation curves; this is the same approach as Hoekstra et al. 
(2001). The second method consists of optimizing the scat- 
ter in the baryonic TuUy-Fisher relation, as in Pfenniger 
& Revaz (2005) and Begum et al. (2008), but interpreted 
in a modified gravity framework (see sect. 14. 2p . We do this 
using the first and the second sample. 

3. Galaxy sample 

3.1. First sample 

We compiled from the literature a galaxy sample of 43 
nearby objects (Table [1]). For each of them, the resolved 
HI kinematical data are available. They are unperturbed 
(without large asymmetries such as lopsidedness or tidal 
interactions) so that a reliable rotation curve can be ob- 
tained. The galaxies have inclinations on the plane of the 
sky that vary from 30° to 90°. A stellar light profile is also 
available in the literature, so that the contribution of the 
visible mass can be computed with reasonable accuracy. 

We have searched for late- type galaxies, where the frac- 
tion of gas (and dark matter) is more important. In partic- 
ular we gathered a large number of dwarf irregulars, which 
are dominated by dark matter. Therefore, they are a con- 
venient laboratory to determine the radial dark matter dis- 
tribution, independently of the stellar mass-to-light ratio. 
FiglH (left) displays a histogram of the morphological type 
in the sample; 70% of the galaxies are late-type. Fig. [2] 
(right) shows the distribution of the maximum velocity of 
galaxies; the median of the sample is Vmax = 128 krn.s"^. 

3.2. Second sample 

A second sample of galaxies was build from the GoldMine 
database. It contains 576 galaxies selected with Vc < 100 
km.s~^. These galaxies are low-mass galaxies with a high 
gas fraction (as in Begum et al,2008). Unlike the first 
sample, they are not kinematically resolved thus the stel- 



lar mass cannot be determined from a rotation curve fit. 
However the GoldMine database provides the stellar lumi- 
nosity in the K-band and the color B-V from which a mass 
to light ratio can be derived (Bell & de Jong, 2001). It 
also gives the width of the HI line obtained by averaging 
the value at 20% of the peak flux with the one at 50% of 
the mean flux and the total mass in HI (see the GoldMine 
website) . 

4. Method 

4.1. Rotation curve analysis 

We analyzed the 43 galaxies (first sample) where a complete 
set of data is available: the luminosity profile, the HI profile 
and the HI kinematics. The atomic gas mass is corrected 
for the primordial helium abundance [Mhc = O.^Mhi). 
The mass-to-light ratios are estimated from the rotation 
curve fits in MOND, by taking into account the stellar and 
atomic disc only (fixing oq = 1.2-10~^° m.s~^). Note that in 
MOND when the baryonic profile is determined the mass- 
to-light ratio is the only parameter of the fit, since the crit- 
ical acceleration oq and the interpolating function /i are the 
same whatever the considered galaxy. These mass-to-light 
ratios are given in Table [21 the stellar mass-to-light ratios, 
in the B-band, vary between 1 and 6 from the late-type to 
early-type galaxies. 

The method consists of modeling all the rotation curves 
to find which pair (oq, c) best fit each galaxy by minimizing 
the . Then we look at the histogram of ao and c to see if 
there is a value common to all the sample. 

If MOND is still correct we expect to obtain a common 
value of ao for all the galaxies. On the other hand, the scale- 
factor c has no reason to be universal. Its value depends on 
the history of each galaxy because of interactions, mergers, 
and the star formation rate. But the mean value should 
be consistent with previous works (Hoekstra et al. 2001, 
Pfenniger & Revaz 2005, Begum et al. 2008) 

4.1.1. Galaxy modeling 

The stellar disc is modelled by an exponential surface den- 
sity: 

i;*(r) = Eo*exp(-r/r*), 

and the bulge is represented by an Hernquist profile, of 
characteristic scale r;,: 



{r/rb){r/rb + 1)^ 

The parameters So^tj t*^ and r^ are deduced from the lumi- 
nosity profile (B-band) fits. 

The HI gas surface density is taken directly from ob- 
servations without any particular modeling. In this way, it 
is particularly constraining to fit the wiggles of rotation 
curves associated with gas overdensity. 

4.1.2. Rotation curve modeling 

The rotation curves of each component are computed us- 
ing the Bessel functions. For a given surface density S, 
the Newtonian circular velocity can be written (Binney & 
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Tremaine, 1994): 

/•OO 

vl{r) = -r S{k)Ji{kr)kdk, 
Jo 

with 

/•oo 

S{k) = -2wG / jQ{kr)T,{r)rdr. 
Jo 

J I and J2 are Besscl functions of order and 1. For an 
exponential disc, it can be simplified to: 

i,l{r) = 47rGSo.r,y2 [I„{y)K,,{y) - h{y)K,{y)] , 

where y = r/2r*. / and K are the modified Bessel functions. 



The MOND rotation curves are deduced from the 
Newtonian acceleration, using the MOND formula: 

aN = aMH{aM/ao) 

and inversely, 

gm = aNv{aN/aQ) 

When not explicitly noted otherwise, we use the standard 
/Lt- function, ijl{x) = x/vT+a?, with x = aAi/ao- 
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4.2. Tully-Fisher relation 



We analyzed the baryonic Tully-Fisher relation using the 
first and second sample, which are independent of that 
of McGaugh et al. (2000) and Begum et al. (2008). The 
method consists of representing cx M with M = M^, + 
cMat and determines which c optimizes the scatter of the 
Tully-Fisher relation. The only difference compared to the 
analysis of the previous authors is that we interpret this 
Tully-Fisher relation in MOND so we fix the slope to rt = 4 
to avoid too many parameters in the fit. 

Contrary to the rotation curve analysis, this one is 
global. This means that we can derive a mean scalefactor c 
between the total gas and atomic gas mass for all the galax- 
ies and not for an individual galaxy. Thus the scatter will 
never be zero. Then when the best value of c for the sample 
is obtained, the critical acceleration is simply given by the 
vertical offset of the relation: log(w'*) — log(M) -I- log(Gao). 

For the first sample the velocity v plotted for the Tully- 
Fisher relation is the last point of the rotation curve. The 
stellar and gas mass is the one computed for the rotation 
curve analysis, integrated until the last point of the rotation 
curve. 

For the second sample, the velocity is obtained from half 
the width of the HI line given in the GoldMine database. 
This velocity is corrected for the inclination (i) of the galaxy 
estimated by the ratio of the short and long axis, 

V = 0.bW/sin{i). 

The stellar mass is computed from the luminosity in the 
K-band. The mass-to-light ratio is given by the color B-V 
(Bell & de Jong, 2001). 



4.3. The visible molecular gas 

In about a dozen galaxies of the sample, there is informa- 
tion on the visible II2 gas, traced by the CO emission. It is 
observed in the more massive galaxies or early-types, due 
to their higher metallicity and consequently greater abun- 
dance of CO. In these galaxies, where the visible molecular 
component is a significant fraction of the gaseous mass, the 
radial distribution of the CO emission is exponential, with 
about the same scale length as the optical exponential disk 
(e.g. Young & Scoville 1991). This radial distribution is 
completely different to the HI distribution, which is more 
extended, with a surface density varying nearly with the 
inverse of the radius. In our approach, we therefore include 
the visible H2 gas in the stellar component; in these mas- 
sive and early-type galaxies, it is never greater than 10% of 
the stellar mass, and enters in the uncertainty on the stellar 
mass to light ratio. 



5. Results 

5.1. The critical acceleration 

In this section, we discuss the dependence of the missing 
mass problem, as a function of scale, on the acceleration. 
Let us reason here in the context of Newtonian gravity. The 
process of rotation curve fitting determines the best mass- 
to-light ratio, thus the visible mass profile of the different 
components (stellar disc, bulge, gas): 

M^,, = Ah + Mfc + M„ 



10^ 



10" 



1 1 

Interpolating function 

\ • 




4-1 


\v«* • - 


■> 

113 
(— 


Newtoi 




9,' 


, .t.rf^---^, 



lO"' 



10-^ 
a/ao 



10" 



10^ 



Fig. 3. In Newtonian gravity, the dynamical to visible mass 
ratio is plotted versus the acceleration of the visible matter. 
At large acceleration, there is no need for dark matter. The 
MOND phenomenology is dictated by the rotation curves, 
which establish a strong correlation between the mass dis- 
crepancy and the acceleration. The gray (color) scale codes 
the type of the galaxies. Late- type galaxies are in dark gray 
(blue) while early- type are in light gray (yellow) . The solid 
line (blue) represents the v function associated with the 
standard /i function (see text). 



and the dynamical mass is defined by 

Mdyn = Mvis + MdM- 

Let consider that the dark matter has a spherical symme- 
try: 

2t \ 2 / N GMoMir) 
r 

When the ratio M^yn/Myis is plotted versus the ac- 
celeration (the Newtonian acceleration of the visible com- 
ponent), a very strong correlation is observed (Fig. [3]). 
Whatever the galaxy considered, the amount of dark mat- 
ter required to fit the rotation curve is the same for a given 
acceleration of the visible component. The color scale indi- 
cates the morphological type of the galaxy. Late-type galax- 
ies are in dark gray (blue) while early-types tend towards 
light gray (yellow). This does not depend much on the es- 
timation of the mass-to- light ratio F^,, as shown by late- 
type galaxies which are gas dominated. This relation has 
also been found by McGaugh (2004), when estimating the 
mass-to-light ratio in several ways (maximum disc, stellar 
populations, ...). Whatever the method, a strong correlation 
between the mass discrepancy and the Newtonian acceler- 
ation is observed. 

This relation can be interpreted in two ways: 

— either we consider the context of Newtonian gravity and 
must find why baryons fall in dark matter haloes so 
that this relation with the acceleration of the visible 
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matter is verified. The ACDM model does not give a 
real clue to this correlation. Van den Bosch & Dalcanton 
(2000) try to obtain this relation using semi-analytical 
models of galaxy formation, in the ACDM framework. 
They must then consider star formation feedback effects 
through supernovae explosions, to reproduce the lack of 
high surface brightness dwarf galaxies. They however 
do not succeed in eliminating massive galaxies with low 
surface brightness. They also tuned their galaxies to fit 
the observed TuUy-Fisher relation in the ACDM model, 
while galaxies automatically reproduce it in the MOND 
model. It should be interesting to look at this in self- 
consistent cosmological simulations with baryons. 
we consider the alternative model of MOND, which is a 
direct application of the mass discrepancy-acceleration 
relation. In this frame, the visible mass is the only one, 
and the gravitation law must be modified to remove the 
need for dark matter, such that Mdyn/ M^s = 1- 

The mass discrepancy-acceleration relation is: 



The mass can be expressed in terms of acceleration: 

GMdyn 



N 
vis 



GMyis ^ / a. 



''vis+DM 



N r I "VIS 

avis J 

ao 



A test particle should feel a MONDian acceleration, a*^ 
equivalent to the Newtonian acceleration with dark matter. 



^vis+DM then. 



a'' = a^isf 



ao 



with the virial theorem, to reveal the mass discrepancy- 
acceleration relation. The TF relation writes as Myis (x V^, 
the virial equilibrium of galaxies as cx Mdyn/ R] from 
the combination of these two relations, it can be derived 
that Mdyn/Myis (X R/V"^ oc ao/a. Although these are only 
global relations, they yield almost the correct slope, ob- 
tained for all galaxies, including points internal to each 
galaxy. The main difference here is that the TF relation 
is an integral part of the MOND model, while it is fine 
tuned in the DM-|-Newton model. 



5.2. Rotation curve analysis 

The minimization of the for the rotation curve fits gives 
the variation of c with respect to oq for each galaxy. Fig. 2] 
(left) displays the mean over the galaxy sample of the factor 
c versus the critical acceleration ap. When ao tends to zero, 
which is the Newtonian limit in the MOND formulation, 
the mean proportion of total gas over the atomic one is 
c — 12±4. In this limit all the dark matter needed to fit the 
rotation curve is in the form of a dark gas disc. This result is 
in agreement with the one found by Hoekstra et al. (2001). 
They argue for c — 1 which is the most common value in 
their frequency histogram but their mean value is larger 
(c = 11, see the error bar in Fig.[51 left). On the other hand 
for oq = 1.2 X 10~^''m.s~^, the actual value of the critical 
acceleration of MOND, the mean c is equal to 1, which 
means that there is no need for additional baryons to fit the 
rotation curve. Thus the critical acceleration ao increases 
as the fraction of dark baryons decreases. This last plot 
gives the general order of magnitude of the values ao and 
c. Now to determine if there is a preferred (universal) value 
of ao, if some dark baryons are taken into account, we plot 
a histogram of ao given by the fit. It is represented in the 
middle of Fig. [4] It can be seen that the critical acceleration 
common to the sample peaks at ao = 0.96 ± 0.39 x 10^^" 



which is equivalent to Milgrom's formulation of MOND: m.s^^. Then, if we look at the histogram of c in Fig [H 

right, it shows that the majority of the galaxy in the sample 
(75%) could contain some dark baryons in the proportion 
1.5 < c < 6, while a few percent of galaxies have a scale- 
factor of 1, equivalent to a model without dark baryons. 
The mean value of the scale-factor is c = 3.2 ± 2.4. Fig. 
[7] displays the different components of the rotation curves 
fitted by MOND with a fraction of dark baryons. 



or 



,N 

ao 



M 



a 
ao 



with 



I(x) — xfi{x) 

v{y) = r\y)/y. 

Thus, rotation curves of galaxies tell us how Newtonian 
gravity should be modified to obtain a ratio Mdyn/ 
of one, without any dark matter halo. The /x-f unction of 
MOND and the critical acceleration ao can be observation- 
ally derived from a Newtonian interpretation of the rotation 
curves. The solid line in Fig. [3] represents the function: 



u{x) = y^o.5 + o.5v/r 



{2lxY 



associated with the standard /i-function. It is in good agree- 
ment with the data. 

Let us note that for any model devised to fit rota- 
tion curves (such as DM and MOND models), it is suf- 
ficient to verify the global TuUy-Fisher relation, together 



5.3. Tully-Fistier relation 

First sample. Fig. [5] (left) shows the evolution of the root 
mean square of the TuUy-Fisher relation fit with respect to 
the fraction between total and atomic gas, c. The factor c 
which optimizes the scatter of the relation is found near 
c = 1.5. This value is lower than for the rotation curve 
method. But this method takes into account only the global 
properties of the galaxies, so it is a statistical study that 
necessitates a large number of candidates. In that case we 
have only 43 galaxies. This is why we analyzed a second 
sample made of several hundred of galaxies. However note 
that if we look at the maximum value of c for which the 
root mean square (rms) of the fit is still better than when 
there are no dark baryons, we find c < 2.5 (see the gray 
region in Fig.[5l left), which is more in accordance with the 
results of the fisrt method. 
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Fig. 4. Left: Variation of the mean c with ag. The hmit oq corresponds to a Newtonian gravity without dark 
matter. The error bar (red) is the result from Hoekstra et ah (2001). When qq 1.2 x 10"^*^ m.s^^, the actual value 
of flo, no additional baryons are required (c 1). Middle: Histogram of aq obtained when fitting the rotation curve 
sample including a dark baryon component. The universal critical acceleration peaks at qq = 0.96 x 10^^*^ m.s^'^. Right: 
Histogram of c when all the rotation curves are fitted fixing qq = 0.96 x 10~^° m.s~^. 



0.32| ^ ^ ^ ^ ^ ^ ^ 1 2.5| ^ ^ ^ ^ 1 2.5 




°'^1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 """-0.4 -0.2 0.0 0.2 0.4 0.6 """-0.4 -0.2 0.0 0.2 0.4 0.6 

c log(v) log(v) 

Fig. 5. First sample. Left: Root mean square of the linear least fit square of the TuUy-Fisher relation as a function of c. 
The scatter of the TuUy-Fisher relation is a minimum for c=1.5. Middle: The baryonic TuUy-Fisher relation for c=1.5 
whatever the galaxy (the masses and velocity are given in the system unit G=l). Right: The baryonic TuUy-Fisher 
relation where the values of c are determined by the rotation curve analysis; they are different from one galaxy to the 
other. 



The baryonic TuUy-Fisher relation where c is fixed to 
1.5, whatever the galaxy, is plotted in Fig[5] (middle). In this 
case, the root mean square is 0.25. To compare this with 
the first method, the same baryonic TuUy-Fisher relation 
is plotted in Fig [5] (right) but the values of c are the ones 
derived from the individual rotation curve analysis, thus 
different in each galaxy. This is why the scatter is lower 
than in the previous plot; the root mean square is 0.14. 

Second sample. With a larger sample of galaxies (n=576) 
selected from the GoldMine database, we performed the 
same study on the baryonic TuUy-Fisher relation. This 
time, the factor that optimizes the scatter of the relation 
is c = 2.6 ± 2.5. The uncertainty corresponds to the range 
of c where the rms is better than if no additional dark 
baryons are considered (Fig.lS]). Moreover, the vertical off- 
set of the relation gives the new value of the critical ac- 
celeration of MOND if some dark baryons are considered : 



ao = 0.85 ±0.35 x 10-^° m.s'^. This analysis of the TuUy- 
Fisher relation is well in agreement with the results derived 
from the rotation curve fits. 

6. Discussion & conclusion 

We revisit the mass discrepancy-acceleration relation, at 
the origin of the motivation for MOND. In the frame of this 
modified gravity theory, this relation is a direct observation, 
which quantifies the critical acceleration and determines the 
interpolation function /i. In the frame of Newtonian gravity 
and dark matter theory, we have to understand the mean- 
ing of this strong correlation, which is related to both the 
TuUy-Fisher and virial relation. Until now, only fine tun- 
ing, involving baryonic physics and supernovae feedback, 
has been invoked as an interpretation. 

We then consider for the first time how the presence of 
dark baryons in galaxies could be made compatible with the 
MOND phenomenology Only a small fraction of all missing 
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Fig. 6. Second sample. Left: Rms of the linear least square fit of the TuUy-Fisher relation as a function of c. The minimum 
is obtained for c—2.6. Right: The baryonic TuUy-Fisher relation including a dark-baryon component, c=2.6 whatever the 
galaxy. 
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Table 2. c is the scale- factor between the total gas (atomic 
and molecular) and the atomic gas. (M^/L)^ is the mass 
to light ratio in the B-band. 



baryons can be present in galaxies, to avoid overpredicting 
the observed rotation curves. If dark baryons are present in 
cosmic filaments in the form of cold gas, it is unavoidable 
to find a small fraction of them in galaxies. We show how 
the presence of these dark baryons in galaxy discs can re- 
duce the critical acceleration uq of MOND. We quantify the 
best pair (oq, c) with c being the scale-factor between the 
total gas and the atomic gas in a galaxy. The main result 
is derived from the rotation curve analysis which involves 
galaxies well resolved with their luminosity profiles and HI 
profiles. More galaxies with high-quality data should be in- 
cluded in the sample, especially the late- type galaxies which 
are very constraining. We applied another method based 
on the minimization of the scatter of the TuUy-Fisher rela- 
tion. We use two different samples, one is the same as for 



the rotation curve fits (with 43 galaxies); the other contains 
a larger number of galaxies (576 candidates). Both meth- 
ods with both samples conclude that the presence of dark 
baryons in the galaxy disc are compatible with MOND. We 
find the scale factor between the total gas and atomic gas 
to be c = 3^2 (by compiling all the results). This factor is 
not expected to be universal because of invidual galaxy his- 
tories. On the other hand, we show that the critical acceler- 
ation of MOND which differentiates the Newtonian regime 
from the mondian regime should be slightly reduced from 
the actual value usually used. In this work the critical ac- 
celeration of MOND, if some dark baryons are taken into 
account, is estimated to beao = 0.96 ± 0.39 x 10^^" m.s~^. 

Our results on the TuUy-Fisher relation are in agree- 
ment with those of Pfenniger & Revaz (2005) and Begum 
et al. (2008), where the error bars are large because of the 
difficulty in determining with very high precision the vis- 
ible baryonic mass. They find respectively c = 3^2 ^'^d 
c = 9^7^- But we interpret the TuUy-Fisher relation con- 
forming to MOND, and we fix its slope to 4. It assumes that 
the velocity used in the TuUy-Fisher relation corresponds to 
the constant circular velocity in the deep MOND regime. It 
is possible that the velocity has not yet reached the plateau 
of the MOND regime. This effect will tend to decrease the 
value of the slope because these samples of galaxies con- 
tain essentially late-types. Their rotation curves are still 
increasing in the outer measured radius, and the velocity 
observed will always be lower than or equal to the asymp- 
totic velocity. 

Let us add that a factor c larger than 7 appears less 
likely from a dynamical point of view (Revaz et al. 2008) . 
With such a large factor, the total mass in gas (HI and 
dark) is about the same order as that in the stellar disc, 
at a redshift z=0. In this case, galaxy discs should be too 
cold and unstable. Note also that the proportion of dark 
baryons estimated in this work corresponds to a quantity 
within ten times the typical scale length of the galaxies, 
due to the extent of the rotation curves. 
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Fig. 7. MOND rotation curve fits: the solid line (red) represents the stellar disc, the atomic gas corresponds to the dash 
line (blue), the dot line (cyan) is the cold and dark molecular gas. The total modelled rotation curve is in dot-dashed 
(magenta) while the observed HI velocity corresponds to the symbols and error-bars. 
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